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Abstract 

We study the low energy behavior of the one dimensional Hubbard model across the Mott 
metal-insulator phase transition in an external magnetic field. In particular we calculate elements 
of the dressed charge matrix at the critical point of the Mott transition for arbitrary Hubbard 
repulsion and magnetization numerically and, in certain limiting cases, analytically. These results 
are combined with a non-perturbative effective field theory approach to reveal how the breaking of 
time reversal symmetry influences the Mott transition. 
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INTRODUCTION 



Exact solutions play an invaluable role in our understanding of the behavior of electron 
systems in one dimension, especially for effective field theory approaches in regimes when 
strong coupling develops. Here effective theories have to be assisted with a non-perturbative, 
unbiased analysis. Strong coupling regimes can develop when the initial interactions between 
electrons are strong, but also in certain situations when the bare couplings are arbitrarily 
weak. One well known example of the latter case is provided by the so called commensurate- 
incommensurate phase transition , in which a single bosonic mode is involved. By 
comparing with a free fermions picture the universal physics at long wavelengths across 
the transition was established A physical system realizing such transition is e.g. the 
repulsive Hubbard model at half filling, where Mott metal-insulator transition takes place 
upon variation of the chemical potential in the absence of an external magnetic field. In this 
case the charge mode undergoes a commensurate-incommensurate phase transition, while 
the spin sector is decoupled. 

In the absence of special symmetries such as time reversal symmetry, however, there is 
no reason for the decoupling of spin and charge modes. This is well known for the Hubbard 
model away of half filling in finite magnetic field s|. At half- filling, on the other hand, 
spin and charge modes are strictly decoupled even in an external magnetic field and the low 
energy sector is equivalent to the Heisenberg antiferromagnet. Here one might be led to 
think that the admixture between spin and charge modes will die out gradually approaching 
the half filling. 

If this expection would be correct, the Mott metal-insulator transition in the presence 
of an external magnetic field would look pretty similar to that at zero field: all relevant 
changes affect only the charge sector, while the spin sector remains practically undisturbed. 
Thus the phase transition would fall in the the single mode commensurate-incommensurate 
universality class even in the presence of magnetic field. In the following we will show that 
the above expectation is misleading. Any small coupling between the spin and charge modes 
increases under the renormalization and qualitatively influences the critical properties of the 
system. 



2 



THE MODEL 



Our starting microscopic model is the Hubbard Hamiltonian of one dimensional lattice 
electrons: 

N N 1 1 

i=l, a i=l 
N , N 

- H^i^i,} + n i,l) ~ g X^-t ~ n i,l) (!) 

i=l i=l 

Here nj )(T is the number operator of electrons with spin a = {f> !} at the i-th site, 4m > is 
the local coupling constant, /i is the chemical potential (half filling corresponds to /i — 0), 
and h is an external magnetic field. 

The asymptotics of the correlation functions of Hubbard model in the presence of two 
gapless modes have been computed using a conjecture for a multivelocity 'conformal' field 
theory 4| by finite size scaling analysis of the low lying energies of ([ID from the Bethe 
ansatz 3j, [5|. Alternatively, the exact finite size spectrum can be used to determine the 
parameters of a mult i- component Luttinger liquid theory in a non-perturbative way. This 
way one obtains an explicit expression for the correct low energy effective field theory for 
the Hubbard model in the sector corresponding to the given filling and magnetic field 6[]. In 
both approaches, the crucial quantity governing the low energy behavior is the k x k dressed 
charge matrix, k being the number of gapless modes in the theory. 

We will first try to follow the second approach and work out the effective field theory 
to better understand non-perturbative effects produced by the magnetic field across Mott 
metal-insulator transition. Then we will compute elements of dressed charge matrix from the 
Bethe ansatz, numerically for generic Hubbard coupling and magnetic fields, and in certain 
limits analytically. Finally we will study response functions and correlation functions and 
discuss the effects of the magnetic field on the Mott metal insulator transition. 

EFFECTIVE FIELD THEORY 
Weak-coupling Bosonization 

First we recapitulate on the effective field theory of the repulsive Hubbard model at half 
filling where umklapp processes open a charge gap and the low energy sector of the Hubbard 



model becomes equivalent to that of the Heisenberg antiferromagnetic chain. When the 
chemical potential is smaller than the gap in the single particle excitation spectrum (charge 
gap) the effective Hamiltonian takes the following form 9|: 

H = | [K s {d x 6 s f + (d x cf> s ) 2 /K s ] + | [K c {d x 9 c f + [d x <p c f/K c ] 

—=fid x (j) c cos (V8iv4> c ) (2) 

'TV 7T Z 



and the only quantities that determine asymptotics of the algebraically decaying correlation 
functions are spin wave velocity v s and Luttinger liquid parameter K s of the spin mode. 
Clearly, spin and charge modes are perfectly separated in (T5]) . However, away of half filling 
and for nonzero magnetic field spin and charge degrees of freedom do not separate any more. 
The reason of the admixture between spin and charge degrees of freedom can be easily 
traced back to the case of noninteracting electrons. For non-zero magnetization, m^O, the 
anisotropy of the Fermi velocities t>f 7^ v± (away of half filling) leads to a coupling of spin 
and charge excitations {(], Q, Q]. For u = this coupling is: 

~ 0t - vi)(d x 4> c d x (j)s + d x e c d x 6 s ). (3) 

For all filling (except of half filling) there is a quadratic coupling between spin and charge 
modes with an amplitude proportional to the velocity anisotropy at u — 0. The modification 
of this admixture for finite u is encoded in the dressed charge matrix. To take fully into 
account this spin-charge coupling one has to work with the renormalized theory with the 
initial free value of the coupling amplitude — v± ~ m(n — l)/2. 



Non-Perturbative Effective Theory 



The effective field theory of the repulsive Hubbard model in a magnetic field can be 
constructed in a non-perturbative way by matching the critical exponents of the multicom- 
ponent Luttinger liquid with those of the Hubbard model 6| in the situation where there 
are two gapless modes. Here, we will carry out this comparison in the limit of the zero 
doping. In the presence of two gapless modes the following asymptotic expression follows 



for the primary fields from the hypothesis of multivelocity conformal theories 



4]: 



(<j>(T,x)<j>(r,x)) 



(v c t + ix) 2Ac (v c r — ix) 2Ac (v s T + ix) 2As (v s T — ix) 2As 



(4) 



where the anomalous exponents are defined in terms of the dressed charge matrix as 5j: 

2A C (A C ) = [z cc D c + Z SC D S + H ZssAN 2 C ~J z CsANs J + 2iV c (iV c ), (5) 
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2A S (A S 



Z CS D C + Z SS D S 



Z CC AN S - Z SC AN C 



' 2 det Z 
At /i = the dressed charge matrix takes the form: 



2N S (N S ). 



(6) 



7 7 
7 7 

^ sc ^ ss 



^/2K C 







(7) 



/K c /2 y/K s /2 

with the familiar Luttinger liquid parameters K c / S of charge/spin sectors. On the other hand 
we have AM = AN C - 2AiV s , 2D C + D s = J c /2k F , D s = -J s /2k F (lOj. Therefore, at zero 
magnetic field one finds A C (A C ) = f(AN c , J c ) and A S (A S ) = /(AM, J s ), i.e. the anomalous 
dimensions of spin and charge fields depend only on spin and charge quantum numbers 
respectively. This is the essence of spin-charge separation [lO]. For finite magnetic field, 
h ^ 0, and below half-filling, n ^ 1, however, this is no longer the case. Here the fields 
that diagonalize the quadratic Hamiltonian (Luttinger liquid fixed point) are not spin and 



charge modes, but rather their linear combinations denoted by 



Heff = v i [(d x <p + ) 2 + [d x e + f] + -± [{d x <p-f + (d x e.) 2 } . 



(8) 



Again, the quantity that connects those fields with microscopic physical spin and charge 
fields is the dressed charge matrix Z: 



detZ 
detZ 



V2 



V2 



(Z s 



+ \Z S c — Z c 



V2 



V2 ' 



(9) 



and 



0_ 



{Z cc Z sc 

{Z s , 



Qc + Qs . ry — Q S 

„ ,e c + e, e c -o, 

Zj rs ) ^= Zl x ~ 



(10) 



V2 ss y/2 ' 

Note that decoupling of the spin and charge degrees of freedom in the effective Hamiltonian 
(IE]) requires <f>± oc <fi c / s and 6± oc 9 C / S , i.e. Z cs = and Z cc = 2Z SC [17|. As we shall see 
below the first condition holds for general values of u and h at half-filling while the second 
is violated for non-zero magnetic field. 



The form (jSJ) of the Hamiltonian is very useful for comparison with the one obtained by 
ordinary weak- coupling bosonization procedure. In particular at zero magnetization where 
spin and charge fields are decoupled we can directly read off Luttinger liquid parameters. 
Indeed, using the explicit expression of dressed charge matrix for zero magnetization in the 
limit of half filling (see matrix (1241) below), the effective theory takes the following form: 



TT _ V C 

Heff ~ ^ 



2(d x <f> c ) 2 + ^(d x 6 c ) 2 



| [(d x <t> s ) 2 + (dJsf] , (11) 



implying that the fixed point values for Luttinger liquid parameters are K c — 1/2, K s = 1. 
Note however, that this theory is valid only at extremely low energies, because the charge 
wave velocity v c — > approaching Mott phase (see Eq. (1291)). 

Once the elements of dressed charge matrix are determined, susceptibilities as well as 
correlation functions can be calculated from the above effective theory. In the following we 
will study the elements of dressed charge matrix at the Mott metal-insulator transition point 
in the presence of magnetic field by Bethe ansatz method. 



DRESSED CHARGE MATRIX AT HALF FILLING 

Within the Bethe ansatz approach the dressed charge matrix is a quantity which deter- 
mines the finite size spectrum and thereby the various susceptibilities for the model. For 
the Hubbard model it is given as: 



(12) 



The elements of Z are solutions of a system of coupled integral equations 
U{k) = 1 + /dA' £c(A')ai(A' - sink) , 

J s 

£ CS (A) = Jdk'co S k'U(k')ai(smk' - A) - ^dA' £ cs (A')a 2 (A' - A) , 
U(k)= /"dA'^A'KtA'-sinAO, 

J s 

£ SS (A) = 1 + Jdk' cosk'U(k')ai(smk' - A) - ^dA' £ ss (A')a 2 (A' - A) . 

Here a n (x) = ^ u^j^ an d j c = J®q, f s = f_ A - In the Bethe ansatz approach the 
boundaries Q and A have to be determined as functions of the filling and magnetization. 

6 



(13) 



At half filling one has Q = ir while A is fixed by the condition ei(±A) = for the dressed 
energy of the magnetic excitations (the chemical potential in the Mott phase is \i — 0): 

ei (A) = h + Au- 4Re v / l - (A - iu) 2 - ( dA' ei(A')a 2 (A' - A) . (14) 

J s 

Using J_ dk cos k /(sin k) = the integral equations (Tl3|) simplify: 

= i , UA) = o , 
Csc(k) = J^dX' £ ss (A>i(A' - sink) , ( 15 ) 

£ SS (A) = 1 - fd\' U(A , )a 2 (A' - A) . 

J s 

With this the dressed charge matrix Z reads at half filling: 

(16) 

'/ i 

sc ss I 

and the remaining entries are determined as 

Z sc = U(k)= fdX' ^(A')ai(A') , 

J. (17) 

Z S s = £ss(^4) • 

For general Hubbard coupling and magnetic field, the integral equations (|14p and ffTo^) 
for ei and have to be solved numerically to obtain the dressed charge matrix ffl6|) as a 
function of the magnetic field or the magnetization. In the weak and strong coupling limits, 
however, one can obtain asymptotic expressions. 

For u — > oo we Jiave the relations Z cc = 1, Z sc = | — ^ for any electron filling n c and 
Therefore we find for strong coupling: 




magnetization m 



a. 



Tfl 

2Z SC = 2 ► 2m for n c -> 1 . (18) 



In the weak coupling regime we analyze the integral equations for half filling (fl4|) and 
f fT5|) for sufficiently small magnetic fields: in this case the boundary A of the remaining 
integration J is large. Rewriting eqs. ( fT5j) as: 

6.(A) = £ + { f° + / A }dA' J R(A-A / )6 s (A / ), 

(19) 



i?(A) 



dcu 



oo 2vr 1 + e 



2u|u| ' 



we can apply Wiener- Hopf techniques to compute £ SS (A) in the relevant region A>i>ii 
where the second integral in the first equation of ffl~9l) can be treated as a perturbation. To 
leading order one obtains: 

1 



y \u) 

G + (uj) 



/ dxe^^ ss (A + x) 
Jo 



-G + (u) 



UOJ 



7TC 



(20) 



Startin g w ith (|20|) higher order corrections to y + (uj) can be computed in an iterative scheme 
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14j giving the diagonal element Z ss of the dressed charge matrix: 



Z ss = lim (-iu)y + (uj) = -= ( 1 

uj—>oo ,/9 V 



Rewriting the expression for Z sc in a similar way we obtain: 



+ . . . 



(21) 



2m 



e, s (A) 



cosh(7rA/2-u) 



1 



7T 



u z — e 2 " + 

2 u \ 2u 1 




(22) 



e 2 



1 

2 " V vre 

Finally, A has to be expressed in terms of the magnetic field or the magnetization. Applying 
the Wiener-Hopf technique to Eq. (Till) the condition e\(A) = gives: 

1 



A 



(23) 



where h c = A^2ir/e Ii(ir/2u) 12|. In particular, we have A — > oo for zero magnetic field 
independent of u. As a consequence the dressed charge matrix at half filling for any u 
becomes: 

1 

1/2 y/2/2 



(24) 



which is the strong coupling limit of the expression obtained in the SU (2) invariant case 
below half filling [5|. For small fields the leading field- dependence of the dressed charge 
matrix is the linear (up to logarithmic corrections): 



2Z SC 



-h. 



(25) 



27t/ 1 (7t/2m) 

Finally, with the known small field behavior of the magnetization m = (h/2irv s ) and the 
spinon velocity v s = 2Ii{iz /2u) / /2u) }12|. this expression can be written as a function 
of the magnetization m: 

2 2tt 



2Z B 



I (tt/2u) 



m 



e -n/2u m for u _^ Q _ 



u 



(26) 




FIG. 1: The combination Z cc — 2Z SC of elements of the dressed charge matrix as a function of the 
magnetization for chemical potential fi = /U_ and coupling as indicated. 

For the crossover between this exponentially suppressed dependence on u to the one 
for the strong coupling (ITS}) the Bethe integral equations (fTI]) and ffT5l) have to be solved 
numerically. The result is shown in Figure [TJ 



SUSCEPTIBILITIES AND CORRELATION FUNCTIONS 

Let us now summarize our results for the elements of the dressed charge matrix for general 
magnetization when doping 5 = 1 — n c — > 0. As we worked out explicitly in this limit: 

(i) Z cc = \/2K c — > 1: this equation simply reflects the fact that holons are extremely 
dilute in n c — > 1 limit and therefore they do not interact directly with each other. 

(ii) Z cs — > 0: this holds because this element of dressed charge is a smooth function of 
density and magnetization, which is well defined at half filling and any magnetization. Due 
to the spin-charge separation it vanishes at half filling. 

(iii) Z ss — > \J K s /2 depends on magnetization and u. In the limit of infinite \u\ it ap- 
proaches the Luttinger liquid parameter of the Heisenberg chain at the same magnetization. 
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For weak coupling \u\ — > we have: 



K s = 1 + — (27) 

TTVf 



where vf = 2 sin ^ is the Fermi velocity of the noninteracting electrons. 

(iv) The last element that we worked out was Z sc . In two limits, namely u — > and 
u — > oo analytic expressions were obtained. The most important result is that, generically, 
Z cc — 2Z SC 7^ in the limit n c — > 1, thus the effective theory does not get completely 
decoupled in spin and charge parts. This has important consequences on response functions 
and correlation functions. It is straightforward to study susceptibilities from our effective 



theory 



151 ]. and indeed one can see that even at the Mott metal-insulator transition point 



the exact expressions for the Hubbard model are recovered: 

I J_ I 4 Zj 

Xc\h=const * Xc\m=const n ry \ 9 \^°J 

71V C 7TV S (1 -2Z sc y 

(these and similar expressions for the spin susceptibilities are given in Ref. [3]). Note that 
the charge susceptibility remains finite across the Mott transition when the magnetization 
is kept constant. 

On the other hand when the magnetization is not kept fixed, then the charge susceptibility 
diverges across the Mott metal-insulator transition, and we can estimate how it diverges with 
the doping. For this we see from Eq. (1281) that it is sufficient to calculate the charge velocity 



as a function of doping. At half filling in the vicinity of the holon band minimum there is 
a relativistic dispersion (Eq. (7.24) from Ref. jl^): E(k) = y/ A 2 + k 2 v 2 (corresponding to 
expansion of the holon band up to second order in momentum) where v is the holon velocity 
at half filling and A is a single particle gap. Therefore, below half filling the charge velocity 
can be obtained asQ]: 



dE 



v 2 8ir , , 

(29) 



k=nS VA 2 + vH 2 8 2 
In particular for u — > oo we have v = 2y/u and it follows v c = 2i\b. 

Once the elements of dressed charge matrix are determined one can easily obtain correla 



tion functions of physical operators 



ll| . We will discuss below only those correlators, 



which decay algebraically on both sides of phase transition. For the Hubbard model those 
correlators involve only spin operators The most slowly decaying one in the Mott in- 
sulator phase is the in plane spin-spin correlation function. In the metallic phase, i.e. the 
phase with two gapless modes, the leading correlator is, depending on the strength of u, 
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either single field (for small values of u) or again in-plain spin correlator (for large values of 
it). The equal time correlation function of in plane spin operator reads: 



G a;a(^>^cr) a .(z| c +Z2J(l/2+l/2(dct Z)2) ( 30 ) 



while in the insulating phase it was given by: 

Ga-A^ < Mcr) = x z c 2 c /2(det Z) 2 = " ( 31 ) 

Here, we can read off the jump in the critical exponent across the Mott metal-insulator 
transition as: 

-=f + *4 + mz? )=1 2 (32) 

Note, that although there is no spin-charge separation the same universal jump of 1/2 in 
critical exponent as in the zero-field case is recovered. 



CONCLUSIONS 



The presence of a magnetic field breaks time reversal invariance and thereby prevents 
the separation of spin and charge degrees of freedom in the ID Hubbard model away from 
half-filling. In this paper we have studied the effect of the resulting admixture of the gapless 
modes on the nature of the Mott metal-insulator transition. We have addressed this problem 
by studying the long wave-length properties of ID Hubbard model across the transition in an 
external magnetic field. Using Bethe ansatz techniques we calculated numerically elements of 
dressed charge matrix for generic cases. In addition exact analytic expressions were obtained 
in the limiting cases for the same quantity. We also constructed a non-perturbative effective 
field theory where the drastic effects produced by the time reversal symmetry breaking across 
the Mott phase transition become manifest: while the susceptibility related to the charge 
mode which becomes gapped at the transition diverges at half filling for fixed magnetic field 
(just as in the time reversal invariant case and as known for the free fermion picture [l| for 
the commensurate-incommensurate transition) it remains finite due to the absence of spin 
charge separation when the magnetization is kept constant. The reason for the behavior 
in the latter case is that there is always an admixture to the effective field of the magnetic 
mode which remains gapless across the transition. 
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One can check that the result (|26p is consistent with a zero magnetization limit of: 
(sin 2irm / y/u cos 7rm)e _7r cos 7rm / 2u 5 obtained after particle-hole transform from the attractive 
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